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Bose-Fermi Mixtures in a Three-dimensional Optical Lattice 
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We have studied mixtures of fermionic ^''K and bosonic ^'^Rb quantum gases in a tliree-dimensional 
optical lattice. We observe that an increasing admixture of the fermionic species diminishes the 
phase coherence of the bosonic atoms as measured by studying both the visibility of the matter 
wave interference pattern and the coherence length of the bosons. Moreover, we find that the 
attractive interactions between bosons and fermions lead to an increase of the boson density in 
the lattice which we measure by studying three-body recombination in the lattice. In our data we 
do not observe three-body loss of the fermionic atoms. An analysis of the thermodynamics of a 
noninteracting Bose-Fermi mixture in the lattice suggests a mechanism for sympathetic cooling of 
the fermions in the lattice. 

PACS numbers: 03.75.Ss, 05.30.Fk, 34.50.-s, 71.10.Fd 
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Quantum liquids and quantum gases are remarkable 
objects vifhich reveal macroscopic quantum phenomena, 
such as superfluidity and Bose-Einstein condensation. 
These fundamental concepts have profoundly influenced 
our understanding of quantum many-body physics. The 
distinct behaviour observed for purely bosonic or purely 
fermionic systems sheds light on the role played by quan- 
tum statistics. New insights can be attained by mixing 
bosonic and fermionic species. One of the most promi- 
nent examples is a mixture of bosonic He and fermionic 
■^He. There it has been observed that with increasing 
admixture of '^He the critical temperature of the transi- 
tion between the superfluid and the normal fluid phase is 
lowered and below the tricritical point phase separation 
is encountered |j]. 

In trapped atomic gases mixing of bosonic and 
fermionic species has led to the observation of interac- 
tion induced losses or collapse phenomena [^ y] and col- 
lisionally induced transport in one-dimensional lattices 
1^. In this work we report on the creation of a novel 
quantum system consisting of a mixture of bosonic and 
fermionic quantum gases trapped in the periodic poten- 
tial of a three-dimensional optical lattice. The optical 
lattice allows us to change the character of the system 
by tuning the depth of the periodic potential. This leads 
to a change of the effective mass and varies the role 
played by atom- atom interactions. The interaction be- 
tween bosonic and fermionic atoms interconnects two sys- 
tems of fundamentally different quantum statistics and 
a wealth of physics becomes accessible which is beyond 
that of the purely bosonic [3, |g or purely fermionic 0, |g| 
case. A variety of theoretical work has been devoted to 
Bose-Fermi mixtures in optical lattices and new quan- 
tum phases have been predicted at zero temperature 
13, llfl lU ll3l • Moreover, the coupling between a fermion 
and a phonon excitation in the Bose condensate mimics 
the physics of polarons 13] . At finite temperature phase 
transitions to a supersolid state and phase separation are 
expected jl4| . 



In our experiment, we prepare fermionic ^''K atoms 
together with a cloud of Bose-Einstein condensed *^Rb 
atoms. The qualitative behaviour when changing the 
mixing ratio between bosons and fermions is depicted in 
figure^ The momentum distribution of the pure bosonic 
sample shows a high contrast interference pattern reflect- 
ing the long-range phase coherence of the system. Adding 
fermionic particles results in the loss of phase coherence 
of the Bose gas, i.e. a diminishing visibility of the inter- 
ference pattern and a reduction of the coherence length. 
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FIG. 1: Interference pattern of bosonic atoms released from a 
three-dimensional optical lattice for varying admixture of Nf 
fermionic atoms at a value Ubb/zJb = 5. The bosonic atom 
numbers are Nb = 1.2 x lO'^ (a and b) and Nb = 8 x 10* 
(c) and the image size is 660/^mx660/im. The coordination 
number of our lattice with simple cubic geometry is 2: = 6. 

Our experimental setup used to produce a degener- 
ate mixture of a Bose and a Fermi gas has been de- 
scribed in detail in previous work [SJ. In brief, fermionic 
^°K atoms are sympathetically cooled by thermal con- 
tact with bosonic ^^Rb atoms, the latter being subjected 
to forced microwave evaporation. The potassium atoms 
are in the hyperfine ground state \F = 9/2, mp = 9/2) 
and the rubidium atoms in the hyperfine ground state 
\F = 2, nip — 2). After reaching quantum degeneracy 
for both species we transfer both clouds into a crossed 
beam optical dipole trap operating at a wavelength of 
826 nm. The laser beams for the optical dipole trap are 
aligned in horizontal plane and their elliptical waists have 
1/e^ radii of approximately 50 /im and 150 /ini in the ver- 
tical (z) and horizontal (x,y) directions, respectively. In 
the optical trap we perform evaporative cooling by low- 



ering the power in each of the laser beams to ~ 35 mW. 
After recompression, the optical dipole trap has the final 
trapping frequencies (w^, , o^y , cj^ ) = 27r x (30, 35, 118) Hz 
for the rubidium atoms. We estimate the condensate 
fraction to be 90% and use this value to obtain the tem- 
perature of both clouds. The Fermi temperature Tp in 
the optical dipole trap is set by the number of potas- 
sium atoms and the trapping frequencies, and we obtain 
T /Tp ~ 0.3, which is in agreement with a direct temper- 
ature measurement of the fermionic cloud. 

The three-dimensional optical lattice is generated by 
three mutually orthogonal laser standing waves at a 
wavelength of A = 1064 nm and a mutual frequency dif- 
ference of several 10 MHz. Each of the standing wave 
fields is focused onto the position of the quantum degen- 
erate gases and the 1/e^ radii of the circular beams along 
the (x,y,z)-directions are (160,180,160) //m. To load the 
atoms into the optical lattice we increase the intensity 
of the lattice laser beams using a smooth spline ramp 
with a duration of 100 ms. This ensures adiabatic load- 
ing of the optical lattice with populations of bosons and 
fermions in the lowest Bloch band only. We have checked 
the reversibility of the loading process into the optical 
lattice by reversing the loading ramp and subsequently 
let the particles equilibrate during 100 ms in the optical 
dipole trap without evaporation. We measure that for 
both the pure Bose gas and the Bose-Fermi mixture the 
condensate fraction decreases by ~ 1.4% per Er of lat- 
tice depth, where Er = h? /2mpibX^ denotes the recoil 
energy and ruRb the mass of the rubidium atoms. 

The physics of the Bose-Fermi mixture in an opti- 
cal lattice can be described by the Bose-Fermi Hubbard 
model (e.g. [lO|). The parameters of the model are the 
tunnelling matrix elements Jb,f for bosons and fermions, 
respectively, and the on-site interaction strength Ubb be- 
tween two bosons and Ubf between bosons and fermions. 
Using the most recent experimental value of the K-Rb s- 
wave scattering length J15| we obtain Ubf /Ubb ~ —2. 

We have studied the phase coherence of the bosonic 
atoms in the optical lattice for various admixtures of 
fermionic particles. We switch off the optical lattice 
quickly and allow for 25 ms of ballistic expansion be- 
fore taking an absorption image of the atomic cloud. 
From the absorption image we measure the visibility of 
the interference pattern. We determine the maximum 
iT-rnax and the minimum rimin of the density of the atoms 
at a momentum \q\ = 2hk with k = 27r/A (see inset 
in figure |3t) [Ul . From this we calculate the visibility 

For the purely bosonic case we obtain results similar 
to previous measurements [la, U^ U^ ■ In our data (see 
figure I^Ji-), the visibility V starts to drop off at a char- 
acteristic value {Ubb/ zJb)c ~ 6.5. For larger values of 
Ubb/zJb the decrease in visibility is approximated by 
V ex {Ubb/zJb)" with v = —1.41(9), which is consis- 
tent with our earlier measurement in a different lattice 
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FIG. 2: a) Visibility of the Bose-Fermi mixture in the optical 
lattice for various mixing ratios between bosons and fermions. 
The intersection between the dashed lines defines the charac- 
teristic value {Ubb /zJb)c- The inset shows the principle of 
the measurement (see text), b) Measurement of the width of 
the central momentum peak which reflects the inverse of the 
coherence length of the gas. The inset shows how the peak 
width is extracted from the column sum of the optical density. 
The dashed line indicates the upper constraint of the width 
imposed by the fitting routine and the error bars reflect the 
fit uncertainty. 



setup giving v — —1.36(5) [la| but different from the 
exponent ly = —0.98(7) obtained in [l2|. For a mixture 
of bosonic and fermionic atoms the results change, see 
also J19j- Whereas in the superfluid regime for very low 
values of Ubb/zJb the visibility is similar to the pure 
bosonic case, the presence of the fermions decreases the 
characteristic value {Ubb/zJb)c beyond which the visi- 
bility drops off significantly. Nevertheless, the visibility 
still shows a power-law dependence on Ubb / zJb with an 
exponent in the range of— l<i^<— 1.5. 

To quantify the shift of the visibility data towards 
smaller values of Ubb/zJb we have fitted the power-law 
decay for large values of Ubb/zJb and extrapolated the 
slope to the visibility for the superfluid situation (dashed 
lines in figure |3^). The intersection defines the charac- 
teristic value {Ubb/zJb)c which depends on the mixing 
ratio between fermions and bosons Np /Nb as shown in 
figure O From this graph it is evident that even very 
small admixtures of fermionic atoms change the coher- 
ence properties of the bosonic cloud significantly. 

The phase coherence of the bosons in the lattice is not 
only described by the visibility of the interference pattern 
but also by the coherence length of the sample 16, 2Q|. 
In a superfluid state the coherence length is comparable 
to the size of the system. The coherence length is related 
to the inverse of the width of the zero-momentum peak 
plus a small contribution from the repulsive interaction 
between the bosonic atoms. For the pure bosonic case 
we flnd that this width starts to increase at a value of 
{Ubb/zJb)c ~ 9 (see Fig. I^b). For the case of a Bose- 



Fermi mixture this value is dramatically altered: for an 
increasing admixture of fermionic atoms to the bosonic 
sample, the value decreases and it behaves very similar 
to the corresponding value (Ubb/zJb)c for the visibility 
(see figure ISl- 
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FIG. 3: Decrease of the coherence of the Bose gas in the 
lattice vs. the admixture of fermions. 



Our interpretation of the simultaneous decrease of 
both the coherence length and the visibility with increas- 
ing admixture of fermions is that the system leaves the 
superfluid phase. While for the pure bosonic case this 
indicates a Mott insulator transition [l^, lul, for the 
mixture the analysis is more delicate due to the differ- 
ent interactions and the different quantum statistics of 
the two species. The full understanding of the observed 
effects including strong interactions and finite tempera- 
ture is challenging. We will consider two limiting situa- 
tions, namely a strongly interacting Bose- Fermi mixture 
at T = in which polarons and composite fermions are 
formed, and a non-interacting mixture at finite tempera- 
ture. In both explanations we encounter a destruction of 
the superfluid with increasing fermionic admixture which 
qualitatively reflects our results. 

At zero temperature several quantum phases of the 
system are predicted [liJ, |lll |lj, |l3, |lj| , depending on 
the sign and the strength of the Bose-Fermi interaction. 
At low depth of the optical lattice the interaction of the 
Bose-Einstein condensate with the Fermi gas leads to the 
depletion of the condensate |21| and to the formation of 
polarons where a fermion couples to a phonon excita- 
tion of the condensate [13 . The couphng strength of the 
fermions to the phonon modes depends on Ubf and the 
ratio Ubb/Jb- If the coupling becomes very strong the 
system is unstable to phase separation {Ubf > 0) or to 
collapse {Ubf < 0). In the stable regime, the polarons 
can form a p-wave superfluid or induce a charge den- 
sity wave, as has been analyzed in one spatial dimension 
|l3J. The enhanced bosonic density around a fermionic 
impurity increases the effective mass of the fermion and 
might enhance the tendency of the bosons to localize. 
For our parameters, the phonon velocity is comparable to 
the Fermi velocity, a regime that is usually unaccessible 
in solids. On the other hand, the interaction of the Bose 



gas with the second species leads to an effective attractive 
interaction between the bosons which would favor a Mott 
insulator transition at a larger depth of the optical lattice 
[lj| . At a larger depth of the optical lattice other effects 
come also into play. Composite fermions consisting of 
one fermion and ub bosons form when the binding en- 
ergy of the composite fermion exceeds the gain in kinetic 
energy that the particles would encounter by delocaliz- 
ing. An effective Hamiltonian for these (spinless) com- 
posite fermions with rcnormalized tunneling and nearest 
neighbor interaction has been derived and their quan- 
tum phases have been investigated theoretically [3, ll3 ■ 
In this situation, the Bose-Einstein condensate can be 
completely depleted by the interactions between bosons 
and fermions. 

For the finite temperature model of the noninteract- 
ing gas we consider the entropy of the cloud of bosons 
and fermions, which is S" = aNpT/Tp + PNB{T/Tcf. 
Tp ~ Tiu)f{&NfY^^ denotes the Fermi temperature 
for Np fermions in a trap with frequency u)f and 
Tc = huJB{NB/C{'^)y^^ the critical temperature for Bose- 
Einstein condensation with a and f3 being numerical con- 
stants. When increasing the depth of the optical lattice 
adiabatically, the temperatures of the two species remain 
equal to each other due to collisions, while Tc and Tp 
evolve very differently. This is due to the fact that the 
tunnelling rates for the fermions are up to an order of 
magnitude larger than for the bosons for our lattice pa- 
rameters. Since the effective masses m|j p oc 1/ Jb,f en- 
ter into the degeneracy temperatures, Tc decreases much 
faster than Tp. At constant entropy, this results in adia- 
batic heating of the bosonic cloud and a reduction of the 
condensate fraction. Simultaneously the fermionic cloud 
is cooled adiabatically, similar to the situation considered 
without a lattice in j22j . For the noninteracting mixture 
with our parameters one expects a reduction of T/Tp by 
a factor of approximately 2 at a lattice depth of 20 Eji- 

In the experiment we have further studied the occu- 
pation of the optical lattice by measuring three-body 
recombination. Lattice sites with a higher occupation 
than two atoms are subject to inelastic losses where a 
deeply bound molecule is formed and ejected from the 
lattice together with an energetic atom. Independent 
of their occupation all lattice sites are furthermore sub- 
ject to loss processes such as off-resonant light scattering, 
background gas collisions, or photo-association due to the 
trapping laser light. The attractive interaction between 
the bosons and the fermions changes the occupation of 
bosons on the sites of the optical lattice. For the given ra- 
tio of the on-site interaction strength of Ubf/Ubb — —2 
it is energetically favorable to have up to five bosons per 
site if a fermion is present. 

The experimental sequence to study the three-body 
decay starts from an initially superfluid Bose gas at a 
potential depth of 10 En. We use a ramp time of 30 ms 
to increase the potential depth of the lattice from zero to 



10 Eji during which we do not observe a loss of atoms. 
Subsequently, we freeze the atom number distribution 
by quickly changing the lattice depth to a large value of 
18 Efi where the tunnelling time of the bosons is tb — 
1/zJb = 23 ms. We monitor the total atom number as 
a function of the hold time in the deep optical lattice 
(see figure 0J and observe two distinct time scales of the 
decay of the atoms. The fast initial time scale is due 
to three-body losses from multiply occupied lattice sites. 
The slower decay is due to single-particle loss processes. 
To extract quantitative information from the loss 
curves we fit the data with the following model. We 
assume that any singly or doubly occupied site decays 
with a single particle loss rate Fi. Multiply occupied 
sites decay with a rate determined by the three-body 
loss constant K^ = 1.8 x 10~^^ cm^/s |2j| and the three- 
body density [n(r)]'^ at the lattice site. Since we start 
from a superfluid the number distribution at the lattice 
sites can be approximated by a coherent state with a 
third order correlation function being equal to unity [2J| . 
We calculate the three-body loss rate assuming gaussian 
ground state wave functions at each lattice site to be 
Fa = 0.24 X n^ s~^, where ns is the number of bosons 
on the site. By fitting the data with this model we extract 
the occupation of the lattice. We obtain ni,2 = 67(3)% of 
the sites with single or double occupation, n^ — 23(9)% 
sites with triple occupation and 714 = 10(8)% of lat- 
tice sites with occupation four. A mean field calcula- 
tion neglecting tunnelling yields the theoretical values 
«i,2 = 58%, fj-a = 33%, and 714 — 17% which gives rea- 
sonable agreement given the simplicity of the model. The 
slow decay rate is determined to be 0.35(7) s~^. 
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FIG. 4: Decay of a pure bosonic gas (squares) and a Bose- 
Fermi mixture (circles) in the optical lattice. The fast initial 
decay of the bosons is much more pronounced in the mixture, 
reflecting the higher density due to Bose-Fermi attraction. 
For the fermions hardly any loss is observed. The error bars 
indicate statistical errors from three repetitive measurements. 

Upon adding fermions to the system we find a much 
faster initial decay due to three-body loss for the rubid- 
ium atoms. The single particle loss constant is, how- 
ever, the same. In contrast, for the fermionic atoms we 
do not observe a particle loss of a comparable order of 
magnitude. This suggests that the observed loss is only 



due to three-body recombination between three rubidium 
atoms. Recent results have suggested that the three-body 
loss constant K^^ for K-Rb-Rb collisions is an order of 
magnitude larger than for Rb-Rb-Rb collisions J3|. This 
is not consistent with our data since we do not observe 
the corresponding fast loss of potassium atoms, similar 
to previous results |25| . 

In conclusion, we have investigated a Bose-Fermi mix- 
ture in a three-dimensional optical lattice. We have ob- 
served that the presence of fermions changes the coher- 
ence properties of the Bose gas and substantially en- 
hances the three-body loss of bosonic atoms. Bose-Fermi 
mixtures in an optical lattice promise to be an incredibly 
rich quantum system [3, llO, U3, llJi llM Q • A number 
of Feshbach resonances between ^^Rb and ''"K [l5|, |2g 
exists which will give access to various quantum many- 
body regimes predicted in the literature as well as to the 
creation of ultracold heteronuclear molecules. 
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